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Abstract: This study aims to simulate the exhaust flow field of 
ship by the method of computational fluid dynamics (CFD) 
concerning with the interference by exhaust temperature, 
shape of stack and rolling angles etc.. In this research wind 
tunnel test for a corvette has been performed to attain 
associated experimental data, which were used as a reference 
basis. During simulation process several configurations of 
stacks have been selected, and combining with various rolling 
angles, exhaust temperatures and velocities, we have gene- 
rated numerous cases from which the diffusion paths and 
temperature distribution of the exhaust flow field can be 
clearly observed and analyzed. In terms of numerical 
simulation, the packaged program computational fluid 
dynamics software has been adopted. The simulation results 
also possess the same trend as the experimental data, which 
have initially confirmed the methods developed here can be 
used for the arrangement of stack and superstructure at the 
stage of initial and conceptual design of ships. 
Keywords — computational fluid dynamics, diffusion 
distribution, wind tunnel test, smoke diffusion, atmospheric 
pressure. 

I. Introduction 

Submarine is the best vehicle to attack surface ships only is 
the biggest threat to maritime lifeline, submarines must 
constantly change shape. When the submarine submarines 
traveling directly affect the shape of the bear resistance, 
especially in the shape of a steep submarine submarines its 
greater resistance line. The streamlined shape of submarines 
(streamlining) submarine -shaped resistance can be reduced, 
however, the factors to be considered in the design also includes 
the cost of loading capacity, regulations and aesthetics. 
Submarine design is hoping to carry more weapons and a 
number of people, so its shape more to the main square, and then 
through the local modification (rounded edge guide, increase 
diversion devices, etc.) in order to reduce drag. How effective is 
this modified, must through analysis of its submarines and thus 
outside the flow and pressure it to calculate the resistance caused 
by the distribution of submarine surface [1,2]. 

The streamlined shape of submarines (streamlining) can 
reduce boat drag, but also need to be considered design factors 
including price, loading, regulations and aesthetics. Such as 
freight wheel design is hoping to be loaded once more and the 
number of goods, so its shape to the main square, and then 
through the local modification (rounded edge guide, increase 



diversion devices, etc.) in order to reduce drag. The effect of this 
revised how must thus calculate the airflow through the analysis 
of its submarines outside its on submarines and surface 
resistance caused by pressure distribution [3, 4]. This ability to 
build a first-class field analysis is an urgent demand for domestic 
shipyards. With the rapid advances in computer technology 
under the new ship launch and shorten the time trend, the use of 
computational fluid dynamics techniques [5, 6] for solving the 
three-dimensional shape of the ship incompressible viscous flow 
field and pressure field in order to establish useful library of 
hydrodynamic and pressure distribution of the modern methods 
used in shipyard design department [7] . 

Value on water dive Motion flow field of calculation 
methods, in the past, mainly in potential theory (potential theory) 
and Reynolds-averaged Navier-Stokes equations governing the 
calculation equations. Potential theory that the central idea of the 
basic assumptions of the ideal fluid flow (ideal flow) and then 
solving the Laplace equation [8], the solution process and the 
establishment of the Green's function integral equation of 
potential flow, this method requires the use of a computer for the 
effectiveness of the more relaxed, so far still widely academic 
use, MC Fang in 2007 this method to solve the catamaran to 
change the unsteady flow field of sports, the results and found 
steady trend in the flow field of high-speed and long-wave 
longitudinal catamaran shake (yaw) and yaw (pitch) has a 
considerable impact [9] ; the same year, Zhu Renchuan, who also 
uses the potential theoretical simulation of fixed offshore oil 
platforms and coastal ship when docked in numerical wave 
flume whole flow field analysis, to be solved with the finite 
volume method [10]; CP Chen in 2005 places the potential 
theory with Dawsonz methods and a small plate method 
catamaran Hing predicted wave resistance of [11], but because 
of a potential theory ignored fluid viscosity of the calculation, 
simulation of fluid motion can be regarded as an idealized 
mathematical model, can't make a better interpretation of the 
fluid flow phenomena for real, the actual occurrence of RANS 
equations on the importance of applying the relative increase in 
many, 2005 FW Hung RANS equations supplemented with VOF 
multiphase flow model for solving the wake characteristics 
submarine movements in unsteady viscous flow field in the 
wake of volatility will increase as the potential depth and speed 
of navigation to reduce the substantial decline [12]; in 2007 N. 
Chang, who use RANS equation k-e turbulence model combined 
submarine jet resolve the situation with the hydrodynamic field 
after the sail boat hull openings, subsequently learned that the jet 



IJER@2015 



Page 118 



International Journal of Engineering Research 
Volume No.4, Issue No.3, pp : 118 - 122 



ISSN:2319-6890)(online),2347-5013(print) 

01 March. 2015 



flow is the main factor affecting the submarine hydrodynamic 
flow field and the field of [13]. 

Computational fluid dynamics calculations applied to the 
hull or underwater vehicle body such as numerical simulation of 
flow field vehicles, generally can be divided into two kinds of 
relative motion with dynamic mesh model; in relative motion 
flow field numerical model, meaning that grid position fixing 
member, leaving the surrounding fluid motion for the case of 
modal flow field analysis, 1996 Y. Tahara, who applied the 
model to solve the boundary ship, tail torrent field and wave 
flow field, experimental places RANS equations with the 
Baldwin-Lomax turbulence model for solving [14]; 2001 CMP 
Sampaio and others to calculate the relative motion model semi- 
displacement boat ilk field of fundamental physical quantity, 
then that is not the same displacement as well as the wedge 
angle on the ship stern body flow field change on very huge 
[15]; Chang Yu in 2006 in unsteady viscous to RANS equation 
solver with the VOF model submarine formation ilk field in 
multiphase flow, the experimental dynamic grid calculation 
domain is divided, waves and surface waves given symmetry 
(symmetrically) boundary conditions at several different Buddha 
labor internally generated and subsequently found to increase the 
number of its internal labor Buddha wavelength also increases 
[16]; the same computing underwater voyage ultra cavitation 
body in motion after the high-speed launch in 2007 YL Xiong 
numerical method and prove such a propulsion unit will be able 
to improve navigation speed at low thrust state [17]; 2002 Ivan 
Malcevic use of dynamic grid computing Lagrangian point in the 
flow, the micro-structure of the blood in which fluid flow field 
desorption [18]. 

In this study, the use of numerical simulation software 
package, will be done with the results of the analysis of the flow 
field measurement results compare to verify the correctness of 
the value of the platform. In this paper, through research, is 
expected to be reached in incompressible flow and mass 
conservation conditions to analyze the flow field containing a 
submerged body skeg for the sail and move them to build "the 
body of the value of the water tank dive" feature. 

II. Governing equations 

Hydrodynamic flow field before carrying value of the 
analysis, we do a few basic assumptions: (1) fluid is a non- 
Newtonian fluid and has a compression (incompressible flow); 
(2) three-dimensional transient flow field flow (transient ); (3) to 
ignore the impact of energy and gravity. 

In order to investigate the hydrodynamic flow field outside 
of submarines, will this paper methods of computational fluid 
dynamics, which means using the finite volume method and 
follow the law of conservation (conservation law) principles for 
the integration of a variety of types that direct discrete physical 
domain its advantages can choose any shape of the grid, 
considering a transfer equation (1) below the form 

M + v.(/^) = v.(rv0) + s, (i) 

dt 

where p the density of water, t the time, V the velocity vector, 

the dependent variable, the diffusion coefficient and the source 
term. Flow field analysis is focused on the solution of three- 
dimensional transient Navier-Stokes equations, which must 



satisfy the continuity equation and the momentum equations 
incompressibility condition of the form are as follows 

8 P 



dt 



+ V-(pV) = 0 



(2) 



dpV 



+ v • (pv ® v ) = -vp * +v • [p eff (vv + (vvy )] (3) 



which jU e g viscosity coefficient ( p e g = p L + p T ) can be 
decomposed into laminar viscosity coefficient /u L , the turbulent 
viscosity coefficient // r for the static pressure P*is defined as 
follows 



P* = P + -pk 
3 



(4) 



For turbulent flow to turbulent eddy viscosity hypothesis defined 
viscosity coefficient ju T . The basic k-e turbulence model 



(5) 



where C M is a constant, k is the turbulent kinetic energy, 8 is the 

rate of turbulent dissipation. Item selection criteria to solve the 
pressure mode, the pressure and velocity coupling term SIMPLE 
algorithm is then chosen, and the choice of second-order line 
item convection upwind difference method, the turbulent flow 
field pattern on the experiment results show that the use of basic 
than the k- 8 model can be more accurate. For k and 8, we can 
separately from the turbulent kinetic energy dissipation rate 
equation and turbulent kinetic energy equation to obtain. 
Turbulent kinetic energy equation: 

+ V • (pVk) = V • [(ju L + ^)Vk] + P k -ps (6) 
Turbulent kinetic energy dissipation rate equation: 

+ V • (pVs) = V • [(p L + i^)Ve] + f (C £l P k - C e2 ps) (7 
dt cj £ k 

) 

where P k = ju T VV • (VV + VV* ) - | V • V(3ju T d -V + pk). 

Within the three-dimensional context, the time -dependent, 
incompressible Navier-Stokes equations are solved together with 
the divergence-free constraint equation using the unstructured 
finite-volume code. In co-located grids, the Rhie-Chow 
interpolation [19] is employed to avoid the oscillatory pressure 
due to even-odd decoupling. For the sake of reducing the 
numerical diffusion error, the second order upwind scheme is 
used to approximate the convection terms, while the diffusion 
terms are discretized by the second-order accurate centered 
scheme. In the present method, a fully implicit backward time- 
stepping scheme for the approximation of temporal derivative 
terms is adopted at the time increment t=0.005s so as to be able 
to predict the true flow physics. To accelerate the convergence, 
the resulting discretized equations are solved within the 
algebraic multigrid context. Finally, some significant physical 
insights will be visualized by virtue of the employed post- 
processor tool. 

III. Verification study 
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Verification of DARPA SUBOFF 

First, in order to prove the availability and accuracy of the 
software packages used in this article, we use a symmetrical 
geometric shape with a DARPA SUBOFF [1] (Figure 1), 
geometry length of 4.356m, the boundary conditions: fluid 
inflow velocity to l.Om/sec, the outlet for the static pressure is 
zero; assuming fluid is water (density 997.0 kg/m 3 , dynamic 
viscosity of 8.899xl0" 4 kg/ms), the Reynolds number Re is 
1.2xl0 6 (in Table 1). Tetrahedron grid system for patterns in 
CFD software automatically generates the unstructured 
computational grid, with a total of 301,971 grid points a. 
Distribution grid as shown, the outer surface of which at 
DARPA SUBOFF, adopt a more dense grid, in order to facilitate 
the analysis of boundary layer flow field. The figure below 
shows the experimental values of the literature [1], numerical 
calculations [2] The results of the comparative analysis. 




(a) 



(b) 



Figure 1: DARPA SUBOFF geometry (a) geometry; (b) surface 
mesh 
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Table 1: DARPA SUBOFF geometry and boundary conditions 

Figure 2 shows the DARPA SUBOFF 0 degrees angle of 
attack pressure C p value distribution, Figure 3 shows the 
numerical difference in the angle of attack for the comparison of 
0 degrees, 8 degrees and 16 degrees when the pressure C p value 
distribution with the literature [2], Figure 4 for DARPA 
SUBOFF respectively, attack angle is 0 degrees, the limit of the 
distribution of surface flow lines 8 degrees and 16 degrees. 
Figure 5 shows the angle of attack respectively 0,8,16 and 30 
degrees to the velocity vector x=2,3 and 4 on the surface of the 
cross section, Figure 6 shows the DARPA SUBOFF attack angle 
of 30 degrees with the surface of the vortex line limits the flow 
line. 

NACA wing flow field verification 

Then NACA0012, 0018 and 0020 respectively monoplane 
placed in numerical hydrodynamic analysis, the geometry shown 
in Figure 7, the wing length of 1.0m, the inflow velocity of 
l.Om/sec, And the computational domain length 10.0 m, 
Reynolds is 7.3653xl0 6 , grid point number is 188,700, 
(tetrahedral grids to 773,898, wedges of 86,150). With 



NACA0012, for example, Figure 8 shows the calculation of the 
surface of the grid NACA 0012 (y =xz plane and 0.5), Figure 9 
for the NACA0012, 0018 and 0020 wing pressure varies with 
the position of the C p value distribution. Figure 10(a) shows the 
NACA0020 Flowing sectional profile, Figure 10(b) is a cross - 
sectional surface of the NACA wing pressure isobars 
distribution [23]. 




:.*.t» .... 4 Mill ua kit 

Figure 2: DARPA SUBOFF 0 degrees angle of attack pressure 
C p value distribution 



CI* Yinger a 




Ottsic* from bcw XI 



(a) 




(b) 



P»t4«nt rasiil 



• Ch ft ige ■■ I 




(c) 



Prwi rrsilt u-=e 





(d) 



Figure 3: DARPA SUBOFF numerical comparison with the 
literature value of the pressure distribution C p 
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x = 30° 














Figure 6: DARPA SUBOFF attack angle of 30 degrees to the 
surface of the vortex line limit flow lines 



Figure 4: Limit the flow line of the distribution of surface, angle 
of attack (a) 16°; (b) 30° 

V. Conclusion 

Flow chart shows also generate a recirculation zone at the 
end of the submarines. Last recirculation zone flow field 
submarine submarines performance impact caused when 
traveling on is great, and it is extremely complex physical 
phenomena, making further understand its mechanism. This 
calculation can capture complex submarine at the vortex 
trajectory. This result can be used in the future to improve the 
appearance of the tail with submarines. Through CFD scientific 
computing (scientific computing), you can understand the 
detailed structure of submarines outflow field distribution, which 
can be used in submarine design process is an important 
reference database. Fully integrated CAD/CAE analysis tools, 
has native build submarines numerical wind tunnel prototype 
platform, through a combination of the future more detail. 




Figure 5: Angle of attack were 0°,8°,16° and 30° degrees 
in x=2,3 and 4 on the surface section of the velocity. 
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Figure 7: NACA wing geometry of the schematic 




Figure 8: NACA0012 calculate the surface of the grid 
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Figure 9: NACA wing of pressure C p value distribution with 
position change 
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Figure 10: NACA 0020 section on the streamline distribution 
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